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Utilizing the Baym-Kadanoff formalism with the polarization function calculated in the random phase 
approximation, the dynamics of the v = quantum Hall state in bilayer graphene is analyzed. Two phases 
with nonzero energy gap, the ferromagnetic and layer asymmetric ones, are found. The phase diagram in the 
plane (Ao,B), where Ao is a top-bottom gates voltage imbalance, is described. It is shown that the energy 
gap scales linearly, AE ~ 14B[T]K, with magnetic field. 



Introduction. — The possibility of inducing and con- 
trolling the energy gap by gates voltage makes bilayer 
graphene [U [21 E] one of the most active research areas 
.with very promising applications in electronic devices. 
Recent experiments in bilayer graphene [H [5] showed the 
generation of gaps in a magnetic field with complete lift- 
ing of the eight-fold degeneracy in the zero energy Lan- 
dau level, which leads to new quantum Hall states with 
filling factors v — 0, ±1,±2,±3. Besides that, in sus- 
pended bilayer graphene, Ref.[3] reports the observation 
,of an extremely large magnetoresistance in the v = 
'state due to the energy gap AE, which scales linearly 
with a magnetic field B, AE ~ 3.5 - 10.5B[T]K, for 
B < 10T. This linear scaling is hard to explain by the 
standard mechanisms JS[ [7] of gap generation used in a 
monolayer graphene, which lead to large gaps of the or- 
der of the Coulomb energy e 2 /l ~ B 1 / 2 , I = {hc/eB) 1 / 2 
is the magnetic length. 

In this Letter, we study the dynamics of clean bilayer 
graphene in a magnetic field, with the emphasis on the 
v = state in the quantum Hall effect (QHE). It will be 
shown that, as in the case of monolayer graphene [8] , the 
dynamics in the QHE in bilayer graphene is described by 
the coexisting quantum Hall ferromagnetism (QHF) [5] 
land magnetic catalysis (MC)[7J order parameters. The 
essence of the dynamics is an effective reduction by two 
units of the spatial dimension in the electron-hole pair- 
ing in the lowest Landau level (LLL) with energy E = 
[9[ [TOl E]. As we discuss below, there is however an 
essential difference between the QHE's in these two sys- 
tems. While the pairing forces in monolayer graphene 
lead to a relativistic-like scaling AE ~ y|eB| for the 
dynamical gap, in bilayer graphene, such a scaling takes 
place only for strong magnetic fields, B > Bthr, where 
our estimate yields B t hr ~ 30 — 60T. For B < B t hr, 
a nonrelativistic-like scaling AE ~ |eB| is realized in 



the bilayer. The origin of this phenomenon is very dif- 
ferent forms of the polarization function in monolayer 
graphene and bilayer one that in turn is determined 
by the different dispersion relations for quasiparticles 
in these two systems. The polarization function is one 
of the major players in the QHE in bilayer, and its 
consideration distinguishes this work from the most of 
previous theoretical ones studying the QHE in bilayer 
graphene [12] 2 \ 

Using the random phase approximation in the anal- 
ysis of the gap equation, we found that the gap in the 
clean bilayer is AE ~ 14B[T]K for the magnetic field 
B < Bthr- The phase diagram in the plane (Aq,B), 
where Ao is a top-bottom gates voltage imbalance, is 
described. These are the central results of this Letter. 

Hamiltonian. — The free part of the effective low en- 
ergy Hamiltonian of bilayer graphene is [T] : 








(1) 



where 7r = p Xl + ip X2 and the canonical momentum 
p = —ihW + eA/c includes the vector potential A 
corresponding to the external magnetic field B. With- 
out magnetic field, this Hamiltonian generates the spec- 
trum E = ±2m, m — ji/2vp, where the Fermi velocity 
vp ~ c/300 and 71 w 0.34 — 0.40cV. The two component 
spinor field &y s carries the valley (V = K, K') and spin 
(s = +, — ) indices. We will use the standard convention: 
^Ks = (ipAi,ipB2)K8 whereas *^, s = (ipB2,ipAi)K's- 
Here A\ and B2 correspond to those sublattices in the 
layers 1 and 2, respectively, which, according to Bernal 
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2 'The polarization effects in bilayer graphene were recently con- 
sidered in 1131 . however, the authors used a polarization function 
with no magnetic field for their estimate. 
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(A2 — B\) stacking, are relevant for the low energy dy- 
namics. The effective Hamiltonian ([IJ is valid for mag- 
netic fields IT < B < B thr . For B < IT, the trig- 
onal warping should be taken into account pQ. For 
B > B t hr, a monolayer like Hamiltonian with linear 
dispersion should be used. 

The Zeeman and Coulomb interactions in bilayer 
graphene are (henceforth we will omit indices V and 
s in the field ^vs)- 

f 9 v q , s e 2 f n o ,n(x)n(x') 

H int = HBBjd 2 xt> + (x)a 3 y(x) + — jd 3 xd 3 x' | X _V 

= [i B B J d 2 x^ + (x)a 3 ^(x) + i y d 2 xd 2 x' [F(x - x') 

x (pi(s)p 1 (s') + p 2 (x)p 2 (x'))+ 2F i2 (x - x')pi{x)p 2 {x')] , 

(2) 

where ps is the Bohr magneton, k is the dielectric con- 
stant, and n(x) = S(z — ^)pi(x) + S(z + ^)p 2 (x) is the 
three dimensional charge density (d ~ 0.3nm is the dis- 
tance between the two layers). The interaction poten- 
tials V(x) and Vi 2 (x) describe the intralayer and in- 
terlayer interactions, respectively. Their Fourier trans- 
forms are V(k) = 2ire 2 /nk and Vi 2 (fc) = 2ire 2 e~ kd / nk. 
The two-dimensional charge densities p\(x) and p 2 [x) 
are: 

pi(x) = * + (x)Pi*(s), ^(z) = * + (x)P 2 *(x), (3) 

where P\ = 1+ | T and P 2 = 1 ~| r are projectors on 
states in the layers 1 and 2, respectively [here t 3 is the 
Pauli matrix acting on layer components, and £ = ±1 
for the valleys K and if', respectively]. 

Symmetries. — The Hamiltonian H = Hq + Hi nt 
describes the dynamics at the neutral point (with no 
doping). Because of the projectors P\ and P 2 in 
charge densities ([3]), the symmetry of the Hamiltonian 
H is essentially lower than the symmetry in mono- 
layer graphene. If the Zeeman term is ignored, it is 
UW(2)s x (2) s x Z$ x Z$, where U^(2) s 

defines the U(2) spin transformations in a fixed valley 
V = K,K', and Z^y describes the valley transforma- 
tion £ — > — £ for a fixed spin s = ± (recall that in 
monolayer graphene the symmetry would be £7(4) |llj). 
The Zeeman interaction lowers this symmetry down to 
G 2 = U { - K \l)+ x x ^"(1)+ x f7( /<0 (l)_ x 

Z 2V x Z 2V' whcre U( v \i) s is the [/(I) transforma- 
tion for fixed values of both valley and spin. Recall that 
the corresponding symmetry in monolayer graphene is 
G x = U { -+\2) v x [/(")(2)y, where U^{2) v is the 17(2) 
valley transformations for a fixed spin. 

Order parameters. — Although the G\ and G 2 sym- 
metries are quite different, it is noticeable that their 



breakdowns can be described by the same QHF and 
MC order parameters. The point is that these G\ and 
G 2 define the same four conserved commuting currents 
whose charge densities (and four corresponding chemi- 
cal potentials) span the QHF order parameters (we use 
the notations of Ref. [5]): 

p s : *J* S = i>lc Aia i>KAis + ^K>A 1 si ) K>Ais 

+ ^KB 2 s^KB 2 s+^K'B 2 s^K'B 2 s, (4) 
p s : *lC* s = i>lc Ma ll>KAi8 ~ i>K> Als 1pK>Ais 

+ ^KB 2 s^KB 2 s-^k, B2S ^K'B 2 s- (5) 

The order parameter (pQ) is the charge density for a fixed 
spin whereas the order parameter ([5]) determines the 
charge-density imbalance between the two valleys. The 
corresponding chemical potentials are p s and p s , respec- 
tively. While the former order parameter preserves the 
G 2 symmetry, the latter completely breaks its discrete 

(s) 

subgroup Z\y. Their MC cousins are 

A s : #+T 3 * s = ^ KAlS i)KA lS ~ 1pK, AlS 1pK'A lS 

- i J KB 2 si J KB 2 s+1pK, B2S i J K>B 2 s, (6) 

A s : $+£t 3 * s = ^ KAi Jka iS + ^ K , AlS i>K>A lS 

~ ^KB^KB^ ~ i>K>B 2 s^K>B 2S ■ (J) 

These order parameters can be rewritten in the form 
of Dirac mass terms [5] corresponding to the masses 
A s and A s , respectively. While the order parameter 
([6]) preserves the G 2 , it is odd under time reversal T 
P3]. On the other hand, the order parameter ([7]) is 
connected with the conventional Dirac mass A. It de- 
termines the charge-density imbalance between the two 
layers pQ. Like p s , this mass term completely breaks 
the Z\y symmetry and is even under T. Note that be- 
cause of the Zeeman interaction, the SU^ v \2)s is ex- 
plicitly broken, leading to a spin gap. This gap could 
be dynamically strongly enhanced |15j . In that case, 
a quasispontaneous breakdown of the SU (i - v \2)s takes 
place. The corresponding ferromagnetic phase is de- 
scribed by p3 = (p + — p~)/2 with the QHF order pa- 
rameter *V 3 *, and by A 3 = (A+ - A_)/2 with the 
MC order parameter ^r 3 a 3 ^ [8]. 

Gap equation. — In the framework of the Baym- 
Kadanoff formalism [TB], and using the polarization 
function calculated in the random phase approximation 
(RPA) , we analyzed the gap equation for the LLL quasi- 
particle propagator with the order parameters intro- 
duced above. Recall that in bilayer graphene, the LLL 
includes both the n = and n = 1 LLs, if the Coulomb 
interaction is ignored [1]. Therefore there are sixteen 
parameters p s (n), A s (n), p s (n), and A s (n), where the 
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index n = 0, 1 corresponds to the n = and n = 1 LLs, 
respectively. The following system of equations was de- 
rived for these parameters: 



(2tt)3 

+ G fsl (a;)k 2 Z 2 /2]y e// (fi-w,|k|) 

e 2 d fi+z . 1-f , 



2nl 2 

G7 s \(n) = sz 1 (n)-i 



2 2 
<i 2 fc 



(8) 



-k 2 i 2 /2 



(2tt) 3 

x k 2 ; 2 /2 + G 5s i(w)(l -k 2 ; 2 /2) 2 
x V e ff (fi-w,|k|) 

e 2 d fi + i , i-e , 



[G{ s0 (w) 



2k/ 2 



(9) 



Here Ai = X)„, s sgn(£L ns ), A 2 = J2 n ,s sgn(£ +ns ), 
and 



1 



UJ — E; 



(10) 

are frequency dependent factors in the bare and full LLL 
propagators, where 

E ins = -(p a (n) + A a (n)) + f(£,(n) - A s (n)) (11) 

are the energies of the LLL states, //o is chemical poten- 
tial, Z is the Zeeman energy, Z ~ /is-B = 0.67-B[T]K. 
The second and third terms on right hand sides of 
Eqs.©, © describe the Fock and Hartree interactions, 
respectively. Note that because for the LLL states only 
the component ipB 2 s (ipAis) of the wave function at the 
K(K') valley is nonzero, their energies depend only on 
the eight independent combinations of the QHF and MC 
parameters shown in Ea.(|lip. The function V e ff(u,k), 
describing the Coulomb interaction, is 



V eff (u,k) 



2ne 2 



1 



-n(w,k 2 ) 



(12) 



where IT(a;,k 2 ) is the polarization function in a mag- 
netic field. Since the dependence of IT(a;,k 2 ) on ui is 
weak, the static polarization will be used. Then, in the 
case of frequency independent order parameters, the in- 
tegration over lo in Eqs. ([§]). ([9]) can be performed ex- 
plicitly, and we get a system of algebraic equations for 
the energies E^ ns of the LLL states. 

It is convenient to rewrite the static polarization 
II(0,k 2 ) in the form II = (m/h 2 )Il(y), where both fl 
and y = k 2 Z 2 /2 are dimensionless. The function Tl(y) 
was expressed in terms of the sum over all the Lan- 
dau levels and was analyzed both analytically and nu- 
merically. At y <C 1, it behaves as H(y) ~ 0.55?/ and 



its derivative IT changes from 0.55 at y = to 0.12 
at y = 1. At large y, it approaches a zero magnetic 
field value, 5(y) — ln4/7r (see Fig.l) 3 ). Because of 
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Fig. 1. The static polarization function 4nll(y). 

the Gaussian factors e - k2( "/ 2 = e ~v m Eqs. © and 
©, the relevant region in the integrals in these equa- 
tions is < y < 1. The crucial point in the analy- 
sis is that the region where the bare Coulomb term k 
in the denominator of V e ff(k) = V e //(0, k) (fT2|) dom- 
inates is very small, < y < 10~ 3 £?[T]. The main 
reason of that is a large mass m of quasiparticles, 
m ~ 10~ 2 m e ~ 10 8 K/c 2 . As a result, the polariza- 
tion function term dominates in V e ff(k) that leads to 
V e ff(k) = C(y)h 2 /ml 2 k 2 , where the part with the fac- 
tor 1/fc 2 corresponds to the Coulomb potential in two 
dimensions, and the function C(y) describes its smooth 
modulations at < y < 1 (see Fig.l). It is unlike the 
case of the monolayer graphene where the effective in- 
teraction is proportional to 1/k. As we discuss below, 
this in turn implies that, in the low energy model de- 
scribed by the Hamiltonian in Eqs. ©, ©, the scal- 
ing AE ~ \eB\ takes place for the dynamical energy 
gap, and not AE ~ y|ei?l taking place in monolayer 
graphene [HI El E] - 

Last but not least, using the model with four- 
component wave functions pQ, we determined the up- 
per limit for the values of B, B t hr-, for which the low 
energy effective model can be used. We found that 
Bthr ~ 30 — 60T, corresponding to the experimental 
values 0.34 — 0.40eV of the parameter 71 = 2mv 2 F . We 
predict that for the values B > Bthr-, the monolayer like 
scaling, AE ~ yJ\eB\ : should take place. 



3 'One can show that the presence of a maximum in the func- 
tion 47rII(y) in Fig. [T]follows from the equality of the polarization 
charge density n(r) in a magnetic field B and that at B = as 
r -»• 0. 
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Solutions. — At the neutral point (/xo = 0, no dop- 
ing), we found two competing solutions of Eqs. ([8]) and 
I) a ferromagnetic (spin splitting) solution, and II) 
a layer asymmetric solution, actively discussed in the 
literature. The energy (fTT]) of the LLL states of the 
solution I equals: 



E 



(J) 



s{Z 



In{B) - 



(13) 



« ns ' 2m/ 2 

where the notation I n (B) is used for the integrals 

Iq(b) = 7^(l + y )e-v = Jdy(l-y + y y-y 

J y/xy + 47rII(y) J ^fxy + 4ttII(#) 

(14) 

with x = 0.003-B(T). Note that the Hartree interaction 
does not contribute to this solution. The situation is 
different for the solution II: 



In{B) 

2ml 2 kI 2 



(15) 



The last term in the parenthesis is the Hartree one. For 
suspended bilayer graphcne, we will take k = 1. 

The energy density of the ground state for these so- 
lutions is (a = It II): 



e<°> = - — V V 



E 



\\Ei a) 



+ (-s0.67B + £A )sgn£g i 



(16) 



It is easy to check that for balanced bilayer (A = 0) 
the solution I is favorite. The main reason of this is the 
presence of the capacitor like Hartree contribution in the 
energy density of the solution II: it makes that solution 
less stable. For Aq = 0, the dependence of the LLL en- 
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Fig. 2. The energies of the LLL states as functions of B. 



ergies E^J S of the solution I on B is shown in Fig. [2] (en- 
ergy gaps are degenerate in £). The perfectly linear form 



of this dependence is evident. Also, the degeneracy be- 
tween the states of the n = LL and those of the n = 1 

LL is removed. The energy gap corresponding to the 

(-0 



v = plateau is AE = (E£>_ - E^ + )/2 ~ 14.3B[T]K. 
In Fig. [31 the phase diagram in the plane (Aq,B) is 




A Q [K] 

Fig. 3. The phase diagram in the (Ao,B) plane. 

presented. The area marked by I (II) is that where the 
solution I (solution II) is favorite. The two dashed lines 
compose the boundary of the region where the two so- 
lutions coexist (the solution I does not exist to the right 
of the dashed line in the region II, while the solution II 
does not exist to the left of the dashed line in the re- 
gion I) . The bold line is the line of the first order phase 
transition. It is noticeable that for any fixed value of 
B (Ao), there are sufficiently large values of Ao (B), at 
which the solution I (solution II) does not exist at all. 
It is because a voltage imbalance (Zeeman term) tends 
to destroy the solution I (solution II). 

In conclusion, the dynamics of bilayer graphene in 
a magnetic field B < B t hr is characterized by a very 
strong screening of the Coulomb interaction that relates 
to the presence of a large mass m in the nonrelativistic- 
like dispersion relation for quasiparticlcs. The func- 
tional dependence of the gap on B in Fig. [5] agrees with 
that obtained very recently in experiments in Ref. [1]. 
The existence of the first order phase transition in the 
plane (Ao,B) is predicted. We also estimate the value 
B t hr, at which the change of the scaling AE ~ \eB\ 
to AE - y/\eB\ occurs, as B thr ~ 30 - 60T. It would 
be interesting to extend this analysis to the case of the 
higher, v = 1,2, and 3, LLL plateaus [H[S]. 
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